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ABSTRACT 

Hidden nonabelian symmetries in nonlinear interactions of radiation with 
matter are clarified. In terms of a nonabelian potential variable, we con- 
struct an effective field theory of self-induced transparency, a phenomenon of 
lossless coherent pulse propagation, in association with Hermitian symmetric 
spaces G/H. Various new properties of self-induced transparency, e.g. soli- 
ton numbers, effective potential energy, gauge symmetry and discrete sym- 
metries, modified pulse area, conserved ?7(l)-charge etc. are addressed and 
elaborated in the nondegenerate two-level case where G/H = SU{2)/U{1). 
Using the f/(l)-charge conservation, a new type of analysis on pulse stability 
is given which agrees with earlier numerical results. 
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1 Introduction 

Since the invention of the laser, much progress has been made in understand- 
ing nonlinear interactions of radiation with matter which made nonlinear op- 
tics a fast developing and independent field of science. Laser light in general 
is expressed in terms of a macroscopic, classical electric field which interacts 
with microscopic, quantum mechanical matter. Unlike the case of classical 
electrodynamics, the electric scalar potential and the magnetic vector po- 
tential do not appear to replace electromagnetic fields in nonlinear optics. 
Instead, the electric field itself, with appropriate restrictions to accomodate 
specific physical systems, plays the role of a fundamental variable which ren- 
ders the problem lacking a field theoretic Lagrangian formulation. However, 
one notable exception is the area theorem introduced by McCall and Hahn 
H]. In 1967, McCall and Hahn have discovered a remarkable coherent mode 
of lossless light pulses which propagate in a resonant, nondegenerate two- 
level atomic medium with inhomogeneous broadening, and named the phe- 
nomenon as self- induced transparency (SIT). In their work, propagation of 
light pulses is depicted in terms of a potential-like variable 6{x), the time area 
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of a suitably chosen electric field, which is determined according to the area 
theorem. Under certain circumstances, the system can be described by an 
effective potential variable ^p{x, t) which satisfies the well-known sine-Gordon 
field theory equation. In this case, the 1-soliton of the sine-Gordon theory is 
identified with the 27r pulse of McCall and Hahn. The cosine potential term 
becomes proportional to the microscopic atomic energy and accounts for 
pulse stability. However, one serious drawback of the sine-Gordon approach 
is its oversimplification. In the sine-Gordon limit, frequency detuning and 
frequency modulation effects are all ignored and microscopic atomic motions 
(inhomogeneous broadening) are not taken into account. Also, the model is 
limited to the nondegenerate two-level case whereas many interesting phys- 
ical systems are degenerate and/or multi- level systems. Until now, only the 
sine-Gordon theory was a known field theory for SIT. Thus, subsequent gen- 
eralizations of SIT to more complex systems have resorted only to the SIT 
equation, i.e. the Maxwell-Bloch equation under slowly varying envelope ap- 
proximation (SVEA), finding soliton type solutions by the inverse scattering 
method in some cases. It should be noted that even though the SIT equation 
expressed in terms of a. U — V pair in the inverse scattering formalism resem- 
bles a nonlinear sigma model equation, it does not provide a field theory like 
the sine-Gordon theory whose equation of motion becomes the SIT equation 
[0]. Following the pioneering work of Lamb Ablowitz, Kaup and Newell 
have extended the inverse scattering formalism to include inhomogeneous 
broadening and obtained exact solutions [^]. In accordance with the area 
theorem, these solutions show that an arbitrary initial pulse with sufficient 
strength decomposes into a finite number of 2ti pulses and Ovr pulses, plus 
radiation which decays exponentially. Extensions to the degenerate as well 
as the multi-level cases have been also found with more complicated soliton 
solutions I, 1,1,0]. 

The purpose of this paper is to clarify nonabelian group structures and 
their physical properties hidden in nonlinear interactions of radiation with 
matter. These group structures allow us to construct an effective field the- 
ory of SIT beyond the sine-Gordon limit. In fact, without inhomogeneous 
broadening all generalizations of the SIT equation which admit the inverse 
scattering formalism can be described by group theoretic effective field the- 
ories which generalize the sine-Gordon theory to the nonabelian cases. We 
extend the potential variable <f to a. matrix valued scalar potential g and 
construct a field theory of SIT in terms of a certain two dimensional non- 
linear sigma model associated with a particular coset G/H for a pair of Lie 
groups G and H. In order to do so, we first find a group theoretic formula- 
tion of the Bloch equation in terms of the motion of a spinning top moving 
along the coadjoint orbit which is determined by G/H, and then express 



3 



the full Maxwell-Bloch equation under SVEA as a field theory generalization 
of the spinning top. We make such a field theory generalization by modi- 
fying the G/H gauged Wess-Zumino-Witten nonlinear sigma model action 
which is a well-known model in mathematical physics. The SU (2)/f/(l) case, 
being identified with the so-called complex sinc-Gordon theory, provides a 
field theory for the nondegenerate two-level SIT. Frequency detuning and fre- 
quency modulation effects are taken into account through the SU{2) valued 
scalar i)()tential g, or stated differently, through three scalar potential func- 
tions. The simplest, sine-Gordon limit arises from the restriction of SU{2) to 
f/(l) X f/(l) and our spinning top interpretation reduces to the well-known 
physical interpretation of the sine-Gordon equation as the equation of mo- 
tion for a chain of pendulums. Other generahzations of SIT to the multi-level 
and the degenerate cases can be made in terms of more complicated coset 
structures. It is shown that all these generalizations correspond to specific 
cosets G/H known as Hermitian symmetric spaces. In the presense of inho- 
mogeneous broadening, we still have the concept of effective potential where 
the matrix scalar potential g is also a function of detuning frequency. That 
is, we have infinitely many microscopic potential functions g labeled by de- 
tuning frequency. This prevents SIT with inhomogeneous broadening from 
a field theory formulation. Certain field theoretic properties, e.g. conserva- 
tion laws, no longer hold in this case which nevertheless receive interesting 
physical interpretations as discussed in Sec. 6. The microscopic potential g 
also reveals the group structure of the dressing method, a solution finding 
technique equvialent to the inverse scattering. 

The effective field theory description of SIT provides a completely new 
angle to the SIT problem and reveals various noble aspects of SIT which 
were not available in the inverse scattering approach. One marked difference 
is the appearance of a group theoretical potential energy term which clarifies 
the topological and the nontopological nature of optical pulses. In particular, 
the periodicity of the potential results in infinitely many degenerate vacua. 
In order for optical pulses to possess finite energy, they have to reach one 
of the degenerate vacua asymptotically as x — > ±00 which assigns specific 
soliton numbers, more than one in some cases, to each pulses. We present 
an explicit form of the potential energy term for various multi-level and de- 
generate cases of SIT and define soliton numbers for each cases. The full 
Lagrangian itself also reveals new types of symmetries of SIT. It possesses 
two kinds of symmetries, first the continuous type; the local H vector gauge 
symmetry and the global U{1) axial vector gauge symmetry, and secondly 
the discrete type; the Krammers-Wannier type duality and the chiral symme- 
try. It is noted that the SIT equation is not invariant under the vector gauge 
transformation and the specific choice of gauge fixing incorporates frequency 
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detuning effect. Tiiis siiows that inhomogeneous broadening is equivalent to 
"averaging" over different gauge fixings. Faraday rotation in tfie presense 
of external magnetic field is also identified with the vector gauge transfor- 
mation. Discrete symmetries generate new solutions from a known one. In 
particular, the Krammcrs-Wannicr type duality relates the "bright" soliton 
with the "dark" soliton of SIT. Our group theoretic formulation of SIT al- 
lows a systematic understanding of the integrability of SIT. We find infinitely 
many conserved local integrals of SIT using the properties of Hermitian sym- 
metric spaces. The group structure of SIT also allows us to find exact soliton 
solutions using the dressing method for SIT. 

The topological nature of solitons is particularly useful in understanding 
the stability of optical pulses. Due to the infinite energy barrier, topological 
solitons are stable against topological number changing fluctuations. In the 
case of nontopological solitons, we show that the U{1) axial charge replaces 
the topological number. By making a pcrturbativc analysis around topologi- 
cal and nontopological solitons and also using the U{1) charge conservation, 
we find the stability behavior of solitons against small fluctuations. This 
result agrees well with earlier numerical work and provides a systematic sta- 
bility analysis based on a field theory formulation which otherwise would not 
have been possible. 

The plan of the paper is the following; in Sec. 2, we introduce a semi- 
classical Maxwell-Bloch description of the optical pulse propagation problem 
and the McCall and Hahn's area theorem. The effective potential concept 
is briefly discussed. In Sec. 3, we construct an effective field theory of SIT 
using a field theory generalization of the spinning top equation. Various 
cases of SIT are associated with specific cosets. Inhomogenous broadening 
is introduced from our field theory point of view and we explain the notion 
of effective potential energy and its degenerate vacua and associated sohton 
numbers. In Sec. 4, the dressing method is explained and applied to obtain 
soliton solutions for the nondegenerate two-level and the degenerate three- 
level system. In Sec. 5, infinitely many conserved local integrals are obtained 
using the properties of Hermitian symmetric spaces. Discrete and continuous 
symmetries of effective theories are aslo explained. The issue of stability is 
addressed in Sec. 6, particularly in regard to the conserved U{1) charge and 
inhomogeneous broadening effect is addressed. Finally, Sec. 7 is a discussion. 
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2 Self-induced transparency 



The multi-mode optical pulses propagating in a resonant medium along the 
a;-axis are described by the electric field of the form, 

m 

E = ^ 'Ei{x, t) exp i{kix — wit) + c.c. (2.1) 
1=1 

where ki and Wi denote the wave number and the frequency of each mode 
and the amplitude vector E/ is in general a complex vector function. The 
governing equation of propagation is the Maxwell equation, 

On the right hand side, electric dipole transitions are treated semiclassically. 
d is the atom's dipole moment operator and the density matrix p satisfies 
the quantum-mechanical optical Bloch equation 

th{^^+v^)p=[{Ho-E.d) , p]. (2.3) 

Hq denotes the Hamiltonian of a free atom and v is the x-component of the 
velocity of the atoms. In most cases, SVEA is further made which assumes 
that the amplitudes E; vary slowly compared to the space and time scales 
determined by ki and wi. Under SVEA, the Maxwell-Bloch equation becomes 
a set of coupled first order partial differential equations for the amplitudes 
E; and the components of the density matrix. Explicit expressions of the 
Maxwell-Bloch equation for some physically relevant cases are given in Sec. 
3. Thus, instead of using quantum electrodynamics for the interaction of 
radiation with matter, a specific choice of physical systems allows an effective 
description using the amplitudes E; as a reduced set of variables. This makes 
the use of potential variables quite difficult if not impossible. Also, lack 
of potential variables causes the physical system to be described only by 
the equation of motion, not by the action principle. Consequently, a field 
theoretic formulation is lacking in pulse propagation problems. However, 
there exists one exceptional case. When pulses propagate in a resonant, 
nondegenerate two-level atomic medium with inhomogeneous broadening, 
McCall and Hahn have introduced an effective potential- like variable and 
shown that an arbitrary pulse evolves into a coherent mode of lossless pulses. 
This phenomenon, known as self-induced transparency, is also observed in 
more general, degenerate and/or multi-level cases. Specifically, McCall and 
Hahn have shown that when the pulse envelope is assumed to be real and 
the time area of the dimensionless envelope function 2E , 

/oo 
dt2E, (2.4) 
-00 
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is an integer multiple of 2tt ( 2nn pulse), then the pulse propagates without 
loss of energy. Otherwise, due to inhomogeneous broadening the pulse quickly 
reshapes to 2mi pulse according to the area theorem, 

— - — = —asmO{x], (2.5) 
dx 

for some constant a. These 2n7r pulses are another type of optical solitons 
which arise from the nonlinear response of matter to the radiation field be- 
sides the most well-known optical solitons of the nonlinear Schrodinger equa- 
tion appearing in optical communication problem. The proof of the area the- 
orem follows from inhomogeneous broadening as well as the Maxwell-Bloch 
equation. The SIT equation for the nondegenerate two-level case is given in 
a dimensionless form by 

dE + 2p <P> = 
dD - E*P - EP* = 

dP + 2i^P + 2ED = (2.6) 

where P is an arbitrary constant and C, = w — wq , d = d/dz , d = d/dz , z = 
t — x/c,z = x/c. The angular bracket signifies an average over the spectrum 
fiO as given by 

<■■■>=/ (■■■)/m. (2.7) 



The dimensionless quantities E, P and D correspond to the electric field, the 
polarization and the population inversion through the relation, 

E = -iE-etod/VQh 

P = -pi2exp[-i{kx-iut)]/AktoNof{0 

D = -{p22-pn)/8ktoNof{0 (2.8) 

where e specifies the linear polarization direction, to is a time constant and 
A^o is related to the stationary populations of the levels.0 In order to under- 
stand the structure of 2mi pulses better, we may impose further restrictions 
such that the system is on resonance (^ = 0) without frequency modulation 
{E being real) and inhomogeneous broadening (/(O = ^(0)- Under such 
circumstances, we can introduce an area function ip{x,t), 

^(x,t) = f Edt (2.9) 



oo 



^For the details of constants, we refer the reader to ref. 
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which in the hmit t ^ oo agrees with 6{x)/2 in (|2.4| ). In terms of ip, the SIT 
equation reduces to the sine-Gordon equation, 



dd(f -2f3sm2ip = 



(2.10) 



together with the identification; 



E = E* = dip , < P >= P = -sin2(/? , < D >= D = cos2<p . (2.11) 
The sine- Gordon equation arises from the Lagrangian 



The periodic cosine potential term exhibits infinitely many degenerate vacua 
and gives rise to soliton solutions which interpolate between different vacua. 
This shows that 2nn pulse is in fact the topological n-soliton solution of the 
sine-Gordon equation in this particular limit. The electric field amplitude E, 
now identified with dcp, receives an interpretation as a topological current. 
Note that the effective potential variable ip is different from the conventional 
scalar or vector potentials of the electromagnetic field. Nevertheless, it is 
remarkable that the potential energy cos2(y9 of the sine-Gordon Lagrangian 
can be identified with the population inversion D which represents the atomic 
energy. Also the Lorentz invariance, which was broken by SVEA, re-emerges 
in the sine-Gordon field theory after the redefinition of coordinates. The 
identification of the atomic energy with the cosine potential term possessing 
degenerate vacua, consequently topological solitons, underlies the stability of 
2mi pulses. 

Unfortunately, the sine-Gordon theory formulation of SIT is too restric- 
tive. In the presense of frequency modulation, E should be complex. There- 
fore, it can not be replaced by a real scalar field p and the sine-Gordon 
limit is no longer valid. Also, frequency modulation effects invalidate the 
area theorem. By a direct application of the inverse scattering, it has been 
found, however, that solitons do exist even in the case of complex E [Q]. 
This suggests that a more general field theory of SIT than the sine-Gordon 
theory could exist in order to account for complex E. Recently, we have 
shown that this is indeed true and the generalizing theory which includes 
both frequency detuning and modulation effects is given by the so-called 
complex sine-Gordon theory in the following way |P; assume that E is com- 
plex and the frequency distribution function of inhomogeneous broadening 
is sharply peaked at ^, i.e. f{C, ) = 6{C, — ^) for some constant ^. Introduce 
parametrizations of E, P and D which generalize ( ^.11| ) in terms of three 
scalar fields ip, 9 and r], 



E = e'^^~^'^\2dr]-^ - idp) , P = ie'^^~^'''Um2p , D = cos2p . (2.13) 




(2.12) 



sm ip 
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These parametrizations consistently change the SIT equation ( p.6|) into a 
couple of second order nonlinear differential equations known as the complex 
sine- Gordon equation; 

— cos 09 — 

ddifi + A — ^T^dridri - 2(3 sm2(p = (2.14) 
sin Lf 

ddii ^ — idndif + dndu)) = (2.15) 

sm 2ip 

and a couple of first order constraint equations, 

2 cos^ Lfdr] — sin^ i^dO — 2^ sin^ = 

2 cos V<9r/ + sin Vc^^ = 0. (2.16) 

Note that the complex sine-Gordon equation reduces to the sine-Gordon 
equation when frequency modulation effect is ignored such that = 0, 9 = 
7r/2 and the system is on resonance (.^ = 0). This reduction is consistent with 
the original equation since solutions of the sine-Gordon equation consists a 
subspace of the whole solution space. The complex sine-Gordon equation 
was first introduced by Lund and Regge in 1976 in order to describe the 
motion of relativistic vortices in a superfluid 0, and also independently by 
Pohlmeyer in a reduction problem of 0(4) nonlinear sigma model [^. The 



integrability and various related properties, both classical and quantum, of 
the complex sine-Gordon equation have been studied since then. In particu- 
lar, the Lagrangian for the complex sine-Gordon equation in terms of (p and 
r] is given by 

S = ^ f dipd^ + 4 cot^ ipd7]dT] - 2(3 cos 2^. (2.17) 
27r J 

This Lagrangian, however, is singular at = nvr for integer n which causes 
difficulties in quantizing the theory. Also, besides the complex sine-Gordon 
equation, the SIT equation comprises the constraint equation ( p.l6|) . Thus 



the Lagrangian ( p.l7| ) does not quite serve for a field theory of two-level 
SIT. In fact, the singular behavior of the Lagrangian ( p.l7| ) is an artifact 
of neglecting the constraint equation. This fact as well as the rationale of 
the above parametrizations can be seen most clearly if we reformulate the 
Lagrangian to include the constraint in the context of a nonlinear sigma 
model as explained in the next section. 



3 Effective field theory 

In order to construct a field theory of SIT in terms of potential variables and 
also find a way to extend to more general multi-level and degenerate cases. 
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we first note that the optical Bloch equation admits an interpretation as a 
spinning top equation hke the corresponding magnetic resonance equations 



np. Denote real and imaginary parts of E and P hy E = E^ + iEi, P = 
Pr + iPi. Then, the Bloch part of the SIT equation ( p.6|) can be expressed 

as 

dS = ClxS (3.1) 

where S = {Pr, Pj, D), Q = {2Ei, - 2En, - 2^), i.e. it describes a 
spinning top where the electric dipole "pseudospin" vector S precesses about 
the "torque" vector Q. This clearly shows that the length of the vector S is 
conserved, 

ISl"^ = P^ + Pf + = 1, (3.2) 

where the length equals unity due to the conservation of probability. The 
remaining Maxwell part of the SIT equation determines strength of the torque 
vector. If Pj = 0, we may solve ( |3.2| ) by taking Pr = — sin 2ip and D = cos 2ip 
and also (p.l|) by taking E = dip as given in (|2.11|) . Then, the Maxwell 



equation becomes the sine-Gordon equation as before. This picture agrees 
with the conventional interpretation of the sine-Gordon theory as describing 
a system of an infinite chain of pendulums. In order to generalize the sine- 
Gordon limit to the complex E and P case, we make a crucial observation 
that the constraint ( |3.2|) in general can be solved in terms of an SU{2) matrix 



potential variable g by 



D P \ _^ /I 



P* -D 
The identity equation. 



g Vg^, ^3 = f Q A ■ (3.3) 



(ap* -od) ^^^^ ^""^^^ ^ ^ ^^'^^ 

gives the Bloch equation if we make an identification, 

g-'dg-R=[^l :|), (3.5) 

where R is an arbitrary matrix commuting with g^^a-^g. The identification 
(|3.3|) requires that Tig^^a^g = and {g^^c^g)^ = (g^^a^g)*. In other words, 
g is a unitary matrix. Similarly, the tracelessness of g~^dg from the identifi- 
cation (|3.5|) further requires g to be an SU{2) matrix. We will show shortly 
that R is determined by demanding a Lagrangian for the above equation and 
is also traceless. Finally, the Maxwell equation becomes 



d{g-'dg - R) 



f -dE 
[dE* 



ip \ . f D P 
-i/s) ' ^ [P* -D 



PWs, g-^o-,g]. (3.6) 
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Thus, we have successfully expressed the SIT equation in terms of potential 
variable g up to an undetermined quantity R. Note that the diagonal part of 
the r.h.s. of ( p.5|) is fixed to a constant accounting for frequency detuning. 
This constrains g such that the diagonal part of g~^dg — Ris equal to (i^, — 
i^). Therefore, for a field theory of SIT, we shall construct a Lagrangian in 



terms of the potential variable g whose equation of motion is (|3.6|) together 
with a Lagrange multiplier for the constraint. In order to help understanding, 
we assume for a moment that R = and the system is on resonance = 0). 



Then, the equation of motion ( p.6| ) arises from a variation of the action 

S = Swzwig) - 5'pot + 'S'const (3-7) 
with the following variational behaviors; 

SgSwzw = J ^'^9{g'^dg)g'^5g, SgS^Q^ ^ 2n I ^"^^^^ ]9'^^9- 

(3.8) 

The action Swzw{g) is the well-known SU{2) Wess-Zumino-Witten func- 
tional, 

Swzw{g) = I Tr g-^dgg'^dg - f Tr g'^dg A g~'^dg A g~'^dg , 

(3.9) 

where the second term on the r.h.s., known as the Wess-Zumino term, is 
defined on a three-dimensional manifold B with boundary S and g is an 
extension of a map (7 : S — > SU{2) to B with g\T, = 9 fill. The potential 
term can be easily written by 

5pot = ^ / "^^a^^g-W (3.10) 

Finally, the constraint requires vanishing of the diagonal part of the matrix 
g~^dg which can be imposed by adding a Lagrange multiplying term >S'(.Qj^g^ 
to the action 



p / TiXa^^g-^dg. (3.11) 
ivr J 



'5'const 2 

The Lagrange multiplier A however adds a new term to the equation of 
motion by 

^s^const = ^ / Tr(-aA(73 + [Aag, g-'dg])g-^5g, (3.12) 

which seems to spoil our construction of a SIT field theory. This problem 
can be resolved beautifully if we lift the action (p.7|) to the "vector gauge 
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invariant" one by replacing the constraint term with a "gauging" part of the 
Wess-Zumino-Witten action, 



5" — Swzw{g) — S-pQ^ + S gauge (3.13) 
^-gauge = ^ J M-Adgg-' + Ag-'dg + AgAg-' - AA) (3.14) 

where the connection fields A, A gauge the anomaly free subgroup U{1) of 
SU{2) generated by the Pauli matrix a^. They introduce a vector gauge 
invariance of the action under the transform 
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h-^gh , A~^h-^Ah + h-^dh , A ^ h'^ Ah + h'^ dh (3.15) 



where h : S U{1). Owing to the absence of kinetic terms. A, A act as 
Lagrange multipliers which result in the constraint equations when integrated 
out. The action ( p.TD may be understood as a particular gauge fixing of the 
vector gauge invariance where A = 0, A = Ka^. However, a more convenient 
gauge fixing which manifests the equivalence with the SIT equation (|2.6| ) is 
where A = i^a^, A = 0. Such a gauge fixing is always possible as a result of 
( |3.23|) . Before proving this, we discuss about the generalization of the action 



( |3.13|) to groups other than SU{2) D U{1). We may simply replace the pair 
SU{2) D U{1) by G D H for any Lie groups G and H and obtain the G/H 
gauged Wess-Zumino-Witten action {Sw zw + Sg^nge) ■ This action is known 
to possess conformal symmetry and has been identified with the action for 
the general G/H coset conformal field theory [|T^. The potential term ( ^.1(J| ) 
breaks conformal symmetry. Nevertheless, it preserves the integrability of 
the model given by (|3.13|) where G/H = SU{2)/U{1) and this model has 



been used in describing integrable perturbation of certain coset conformal 
field theories fl^, [1^. For a general pair of G and if, the expression for 



the potential which preserves integrability has been also found |T^, TB|. It is 
given by 



Vt = '^'^aTg-'T (3.16) 

where T and T are constant matrices which commute with the subgroup 
if, i.e. [T, }i\ = [T, h] = 0, for h E H so that the potential term is vector 
gauge invariant. In general, SpQ^^ is specified algebraically by a triplet of Lie 
groups F D G D H for every symmetric space F/G, where the Lie algebra 
decomposition f = g © k satisfies the commutation relations, 

[g , g] C g , [g , k] C k , [k , k] C g . (3.17) 

We take T and T as elements of k and define h as the simultaneous central- 
izer of T and T, i.e. h = Cg(T,f ) = {B e g : [5 , T] = = [5 , f]} 
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with H its associated Lie group. With these specifications, the action ( p.l3| ) 
becomes integrable and generahzes the sine-Gordon model according to each 
symmetric spaces [jl^. For compact symmetric spaces of type II, e.g. sym- 
metric spaces of the form G x G/G, the elements g and T take the form 
g ® g and T ® 1 — 1 ® T (and similarly for T). In which case, the model be- 
comes effectively equivalent to the case where T, T belong to the Lie algebra 
g. Thus the model is specified by the coset G/H where H is the stability 
subgroup of T and T for T, T G g. In this paper, we will restrict ourselves 
only to this case. As we will see later, physically interesting cases of SIT 
all correspond to even more specific symmetric spaces, where G/H becomes 
Hermitian symmetric spaces and the adjoint action of T defines a complex 
structure on G/H. 

Now, we demonstrate the integrability of the model by expressing the 
equation of motion arising from the action ( |3.13| ) in a zero curvature form, 

5gS = hiiXd + g-^dg + g-^Ag, d + A]+P[T, g-'fg])g-^5g 

Zn J 

= fTr[d + g-'dg + g-'Ag + p\T , d + A+ \g-'fg ]g-'Sg 

Zn J A 

= (3.18) 

where A is an arbitrary complex constant. That is, the equation of motion is 
given by a zero curvature condition in terms of a U — V pair, 

[d-U, d-V]=0, (3.19) 

where 

U = -g~'dg - g~'Ag - p\T, V = -A- \9-'fg. (3.20) 

This shows that the equation of motion arises as the integrability of the 
overdetermined linear equations; 

{d + g-'dg + g-'Ag + /5AT)M/ = 0, {B + A + jg-'fg)^ = 0. (3.21) 

The constraint equations coming from the A, A- variations are 

SaS = ^jTTi-dgg-' + gAg-'-A)5A = 

6aS = ^lTTig-'dg + g-'Ag-A)6A = Q. (3.22) 

Note that these constraint equations when combined with ( p.l^ ) imply the 
flatness of the connection A and A, i.e. 

F,s=[d + A, d + A] = . (3.23) 
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In the following, we show that ( p.l8 ) and ( 3.22|) can be identified with the 



SIT equation for various cases depending on the choices of the groups G and 
and T, T as well as the specific choice of gauge fixing. Thus, we obtain 
the action principle of SIT equations when inhomogeneous broadening is ig- 
nored. The inclusion of inhomogeneous broadening is also obtained in our 
field theoretic context. Remarkably, in the presense of inhomogeneous broad- 
ening, the notion of effective potential still persists and the inhomogeneous 
broadening effect, i.e. Doppler shifted atomic motions, can be beautifully 
incorporated into the U{1) vector gauge transformation of the theory. 

3.1 Examples 



Nondegenerate two-level system 

This is the simplest case of SIT which was originally considered by McCall 
and Hahn. It also accounts for the transitions 1/2 1/2, 1^0, 1^1 and 
3/2 ^ 1/2 for linearly polarized waves and the transitions 1/2^1/2, 1^0 
and 1^1 for circularly polarized waves. The SIT equation is given by (|27 
which can be expressed in an equivalent zero cuvature form. 



i^X + ii -E \ r._l_(D P 

E* \\P* -D 



. (3.24) 



In order to show that this SIT equation in fact arises from the effective field 
theory (CT) , we take H = U{1) C SU{2) = G and T = -f = ia^ = 
diag(2, —i) for Pauli matrices (Tj. We fix the vector gauge invariance by 
choosing 

A = i^a-i , A = (3.25) 

for a constant Such a gauge fixing is possible due to the fiatness of A, A. 
Comparing ( |3.20| ) with ( |3.24| ), we could identify E,P and D in terms of g 
such that 

_i / -E\ .{DP 



9''dg + ig-'Tg-iT=[ 1 , g-'Tg = -i{p, (3.26) 



which are consistent with the constraint equation (|3.22|) . Comparing with 



(|3.3|)-( |375D , we find that R is fixed to —iE,g ^Csg. Then, the zero curvature 
equation ( 3.18|) agrees precisely with (|3.24|) . Furthermore, if we make an 



explicit parametrization of the 2x2 SU (2) matrix g by 

tsmipe~'" e~""'cos(y9 



„ _ gjr?(T3gi(^(cos 6*0-1 -sin 6*0-2 )gjr?(T3 — ( ^ COS (fi 2 Sm (^6 \ f3 27) 
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we recover the parametrizations of E,P and D as given in ( p. 13 ) and the 



SIT equation becomes the complex sine-Gordon equation ( ^.141 ) ( PTT^D and 



the constraint equation ( p.l6| ). The potential term in (|3.13| ) now changes 
into the population inversion D, 

Spot = /^cos2^ = J(^D, (3.28) 

which for /? > possesses degenerate vacua at 

f = fn = + 2^'^^ n & Z and 6 = 6q for 6q constant . (3.29) 

The property of degenerate vacua and the corresponding soliton solutions 
will be considered in Sec. 4. 

Degenerate two-level system 

One of the deficiencies of the SIT model of McCall and Hahn is the 
absence of level degeneracy. Since most atomic systems possess level degen- 
eracy, the analysis of the nondegenerate two-level system does not apply to a 
more practical system. Moreover, level degeneracy in general breaks the inte- 
grability and does not allow exact soliton configurations. For example, prop- 
agation of pulses in a two-level medium with the transition jh = 2 —>■ ja = 2 
is effectively described by the double sine-Gordon equation 

ddip = Ci sin (p + C2 sin 2{p (3.30) 

which is not integrable. Nevertheless, there are a few exceptional cases which 
are completely integrable even in the presense of level degeneracy. It was 
shown that [§, ^ the SIT equations for the transitions jb = ja = 1, jb = 
1 ja = and jfc = 1 — s> ja = 1 are integrable in the sense that the SIT 
equations can be expressed in an inverse scattering form. In the following, 
we show that these cases correspond to the effective theory with G = SU (3) 
and H = U{2) <Z G. Also, we show that the local vector gauge structure 
incorporates naturally the effects of frequency detuning and longitudinally 
applied magnetic field. Consider a monochromatic pulse propagating through 
a medium of degenerate two-level atoms in the presence of a longitudinal 
magnetic field. Then, the Maxwell-Bloch equation under SVEA is given by 

fj,m 



1 \ 
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(3.31) 



qm 



The dimensionless quantities e'' and R are propotional to the electric field 
amplitude E and the density matrix p, where q is the polarization index 
and the subscripts fi, fi , . . . and m, m , . . . denotes projections of the angu- 
lar momentum on the quantization axis in two-level states \a > and \b > 
respectively.0 J denotes the Wigner's 3j symbols 



\Ob-rn 



V3 



3a 1 3b 

—m q 



(3.32) 



and VLa{VLh) is a dimensionless coupling constant of an external magnetic field. 

In general, ( p.31| ) is not integrable. However, with particular choices ja 
and jbi ( |3.31| ) can be recasted into the zero curvature form, or the U ~V 
pair as in (|3.19|) . Specifically, for the transition = 1/2 ^ = 1/2, U and 
V becomes 



U 



U+ 
f/_ 



V 



V+ 

v_ 



(3.33) 



where 



~i{x + A) ±ie 



±1 



±1* 



i{x + A) 



V+ = —. 



1 -r- 1 1 



+ 2^2 
(a) 



1 



+ 2^2 



^2^2 ■ 



X 



(3.34) 



In the context of effective field theory, we again identify the U — V pair in 
terms of g by 



U = -g-'dg - g-'Ag - f3\T, V = -jQ-'Tg 



where the gauge choice is 



A 



and 



/ —ix 

ix 



V 

T = -f 



0\ 



-ix 

ix J 



(73 

as 



A = 0, 



(3.35) 



(3.36) 



(3.37) 



^For details of propotionality constants and their physical meanings, we refer the reader 
toref. B. 
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with the Pauh matrix a-^. Here we set (3 = 1 for convenience. The resulting ef- 
fective theory is specified by the coset G/H = iSU{2) x SU{2))/{U{1) x f/(l)) 

with gi, g2 C SU{2) and the two U{1) subgroups 





such that g ■ 
are generated by 



g2 
as 




and 



Note that the specific form of 



^3. 

the identification ( p. 341) requires gi and g2 to be SU (2) as in the case of the 
nondegenerate two-level system. Thus, this case is simply two sets of the 
nondegenerate two-level case. 

Another integrable cases are the transitions; = 1 ^ ja = and jb = 
^ ja = 1. In each cases, the U — V pair is given by 



U 



X 



X 



and 



V 



2A 



V = — 



2A 



/ -^00 

p(M* 



A + i{x - 




—is 


-1 








f.A- 


ix 













I' 


y 






for 


jb 


y 






for 


jb 




■ 1 

0. 



(3.38) 



-fX_l_l 
p(a) 

jj(ba) 

pW 

-fXl_l 



p(M 

pW 



for jb 



for jb = l^ Ja = 0. (3.39) 



-1-1 



The gauge fixing is given by 



A 




A = 0, 



(3.40) 



and the effective field theory is specified by G/H = SU(3)/U{2) with 

/2 



T = -T 



2i 
3" 





Vo 



-1 



(3.41) 







and (3 = 1. It is interesting to observe that the electric field components 



e ^ in ( p.38[ ) parametrize the coset SU (3)/f/ (2) and the vector {e e^*Y 
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transforms as a vector under the U{2) action. In particular, since frequency 
detuning amounts to the global f/(l)(c f/(2)) action while longitudinal mag- 
netic field amounts to the global f/(l) x t/(l)(c U{2)) action, the effects of 
detuning and magnetic field to e'^,e~^ can be easily obtained. 

Three level system 

The propagation of pulses in a multi-level medium with several carrier 
frequencies as given in ( |2.1D is a more complex problem than the two-level 
SIT and in general is not an integrable system. However, with certain re- 
strictions on the parameters of the medium, it becomes integrable again and 
reveals much richer structures. Here, we will not exhaust all integrable multi- 
levels cases but restrict only to the degenerate three-level case and provide 
an effective field theory formulation. Other multi-level cases can be treated 
in a similar way. The propagation of two-frequency light interacting with 
a three-level medium, having either the A or V^-type resonance configura- 
tions, has been studied earlier and its integrability has been demonstrated 
in the context of the inverse scattering method ||^. We suppress the gen- 
eral Maxwell-Bloch equation formulation for the three-level case and refer 
the reader to Ref. for details. Here, we extend the Maxwell-Bloch equa- 
tions of Ref. [0 in order to include a longitudinal magnetic field. Then, the 
Maxwell-Bloch equation in a dimensionless form, describing the A configu- 
ration with jb = 1, ja = jc = 0, is given by 

%^ = -zp^^ J = 1,2, q = ±l (3.42) 

and 

id + ttoihv-2A,-n,q))pl = -t(J2eim^,^-eln,-elr) 

I 

Q 

id + itoik2V-2A2-n,q))pl = -li^eim.^-eln^-eY) 

I 

(a-zto(M-A;it;-2A2 + 2Ai))r = -2 ^(^^^2* - 

<? 

{d + im,{q-q))m^^, = ^ (^f - e/pf ) (3.43) 

i=i,2 

where e'j, j = 1, 2 is the amplitude of a double- frequency ultrashort pulse 
and q = ±1 denote the right (left) -handed polarization. Other variables are 
proportional to the components of the density matrix 

Pi = P-qo exp[-i{kix - wit)]/Na, pI = exp[-i(A;2X - W2t)]/Na 
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r = -PocT^ exp[i(/i;i - k2)x - i{wi - W2)t]/Na (3.44) 



and to is a constant with the dimension of time and Na is the polulation 
density of the level \a >. 2Ai = wi — Wba, 2A2 = W2 — Whc measure the 
amount of detuning from the resonance frequencies. The integrability of 
( |3.43|) comes from its equivalent zero curvature formulation with the 4x4 
matrix U — V pair, 

^ ^ /-Ai-iAl2x2 -iE 

V -iE^ -A2 + i\l2x2 

^ = ^ ( Pt _n) ^3.45) 

where 



\ e{ 62 J \ Pi P2 J V "^-11 iTT'ii 

with 

a = ^-^{kiv + k2V - 2Ai - 2A2 + AQb) 

b = ^(fcif + M-2A1 -2A2-4fib) 

X = ^-^(-Skiv + k2V + QAi - 2A2) 

y = '-jikiv - 3k2V - 2Ai + QA2) (3.46) 

In the context of effective field theory, this corresponds to the case where 
G/H = SU{A)/S{U{2) X U{2)) with the gauge fixing 

A=(j' ^,),A = (3.47) 



and 



A2 



T = -T = t{^'^' _^ ). (3.48) 



2x2 



Similarly, we may repeat an identification for the case [|^, ja = jc = 1, jfe = 0, 
and can easily verify that it corresponds to the symmetric space SU (5)/?7(4). 
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Thus, we have shown that various cases of SIT all correspond to the effective 
field theories with specific symmetric spaces, i.e. 



SU{2)/U{1) ^ nondegenerate two-level 
SU{3)/U{2) ^ degenerate two-level; j;, = ^ ja = 1, 

= 1 3a =0, jb= l-^ 3a = 1 

(Sf/(2)/f/(l))2 ^ degenerate two-level; = 1/2 ^ ja = 1/2 
SU{A)/S{U{2) X U{2)) ^ degenerate three level; ja = 3c = 0, jb = 1- 
SU{5)/U{4:) ^ degenerate three level; ja = jc = 1, jt = 0. 

(3.49) 

Each case corresponds to a special type of symmetric spaces known as Her- 
mitian symmetric spaces. Hermitian symmetric space is a symmetric space 
equipped with a complex structure which in our case is given by the adjoint 
action of T up to a scaling. In Sec. 5, the characteristic properties of Hermi- 
tian symmetric space will be used in obtaining infinitely many conserved local 
integrals. Our association of various SIT systems with Hermitian symmetric 
spaces suggests that to each Hermitian symmetric space there may exist a 
specific SIT with a proper adjusting of physical parameters. Especially, the 
multi-frequency generalization of SIT in a configuration of the "bouquet" 
type may correspond to the symmetric space SU{n)/U{n — 1) for an integer 
n. However, for large n, it requires a fine tuning of physical parameters which 
makes the theory unrealistic. 

3.2 Inhomogeneous broadening 

So far, we have identified various cases of SIT with specific effective field 
theories when inhomogeneous broadening is neglected. In the presense of 
inhomogeneous broadening, due to the microscopic motion of atoms, each 
atom in a resonant medium responds to the macroscopic incoming light with 
different Doppler shifts of transition frequencies. Thus, microscopic variables, 
e.g. the polarization P and the population inversion D are characterized by 
Doppler shifts and they couple to the macroscopic variable E through an 
average over the frequency spectrum as given in ( p.7|) . A remarkable property 
of our effective field theory formulation is that it incorporates inhomogeneous 
broadening naturally only with minor modifications. The notion of effective 
potential variable g still persists where the microscopic variable g becomes 
also a function of frequency ^, i.e. g = g{z,z,C,)- However, the action 



principle in ( |3.13| ) is no longer valid despite the use of the potential variable 
g. In order to formulate the SIT equation with inhomogeneous broadening 
in terms of the potential variable g, we relax the constraint equation ( p.22| ) 
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and require only 

{g-^dg + g-^Ag)^-A = Q, (3.50) 

where the subscript specifies the projection to the subalgebra h. The hnear 
equation is given by 




d + g-^dg + g-^Ag -^T + XT)^ = 

g-'fg' 



^ = (3.51) 



where the constant A is a modified spectral parameter which becomes A + ^ 
in the absense of inhomogeneous broadening. The angular brackets denote 
an average over ^ as in ( ^.7| ). As in the previous examples without in- 
homogeneous broadening, we make the same identification of the matrix 
g~^dg + g~^Ag — with various components of macroscopic electric fields 
which are independent of the microscopic quantity ^. This requires the ^- 
dependence of g{z, z,^) to be such that g~^dg + g^^Ag — S,T is independent 
of ^. It is easy to see that this requirement is indeed satisfied by various inte- 
grable SIT systems we have considered previously except for the three-level 
system, where we must take ^ = — to^i = — ^o^2- It means that in order 
to preserve the integrability in the presence of inhomogenuous broadening, 
two detuning parameters of the three-level system must be equal. Note that 
\E'(A, z) is not a function of ^. The integrability of the linear equation ( p.51j ) 



becomes 







d + g-'dg + g-'Ag-iT + \T , 9 + (t^ 
-d{g~'dg + g-'Ag-^T) + ([T, g-'fg]) (3.52) 



where we used the fact that g ^dg + g ^Ag — is independent of ^ and also 
the identity 

d{g-'fg) + [g-'dg + g'^Ag, g^'fg] = 0. (3.53) 

Identifying g~^Tg with components of the density matrix as in the previous 
cases, we obtain the SIT equation with inhomogenous broadening. For ex- 
ample, we may identify E, P and D as in ( p.26| ) so that ( |3.52|) and ( p.53| ) 



become the SIT equation with inhomogeneous broadening for the nondegen- 
erate two- level case given by (p.6|) . Note that each frequency ^ corresponds 
to a specific gauge choice of the vector U{1) subgroup. Therefore, inhomoge- 
nous broadening is equivalent to averaging over different gauge fixings of 
f/(l) C H. This implies that inhomogenously broadening can not be treated 
by a single field theory thereby lacking a Lagrangian formulation. It is re- 
markable that nevertheless the group theoretical parametrization of various 
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physical variables in terms of an effective potential g is still valid. This group 
theoretical structure underlies the integrability of the SIT equation as we will 
show in Sec. 4. 



3.3 Effective potentials 

One of the main advantage of our effective field theory approach to SIT 
problems is the introduction of a notion of effective potential energy. This 
allows us a systematic understanding of stability of solutions as well as the 
topological nature of soliton solutions. The effective potential energy term 
in (|3.16|) reveals a rich structure of the vacuum of the theory. In general, the 
potential (|3.16| ) is "periodic" according to the coset structure G/H. This 
results in infinitely many degenerate vacua, which are specified by a set of 
integers, and also soliton solutions - finite energy solutions which interpolate 
between different vacua. In the nondegenerate two-level case, the potential 
term ( |3.1CI| ) becomes a periodic cosine potential ( ^.281 ) and each degenerate 
vacua are labeled by an integer n as in ( |3.29| ). Solitons interpolating between 
two different vacua with labels Ua and nt as x varies from — oo to oo are 
characterized by a soliton number An = nb — Ua. In order to understand the 
vacuum structure of the potential for a more general G/H case, we first note 
that the potential term Ti{gT g~^T) is invariant under the change g —>■ gh for 
h E H. Consequently, we may express the potential term by a coset element 
m G G/H such as Tr(mTm~^T), where m is given specifically by 



/ B\ ( cos/Bfit B^WB \mVBW' 



m = exp 



V-5t Oj {-sinVB^VWB ^B^ cos^/WB 

(3.54) 

The matrix B parametrizes the tangent space of G/H. This manifests the 
periodicity of the potential through the cosine and the sine functions. In the 
case of our interest, i? is a complex- valued matrix of dimension 1x1, 1x2 and 
2 X 2 for the coset SU{2)/U{1), SU{?>)/U{2) and SU{A)/S{U{2) x U{2)) re- 
spectively. Then, due to the relation: B sin V B^B^/ B'^B = sin V BB'^V BB^ 
the potential term reduces to the form 



Tr ( / — 2 sin VBB^) + Tr (/ - 2 sin^ /bTb) (3.55) 

for SU{2)/U{1) and SU{A)/S{U{2) x U{2)) cases and 

Tr (4/ - 6 sin^ V^fit) + Tr (/ - 3 sin^ Vst^) (3.56) 

for SU{3)/U{2) case. In order to reduce further, we define non-zero eigen- 
values of B^B by = l,..,r = rank{i?^i?}, which are positive definite 
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and coincide with those of BBK With these eigenvalues, the potential term 
takes a particularly simple form, 

-6^sin2 0,, (3.57) 



a 



where the positive constants a and b can be read directly from ( |3.55D and 
( |3.56| ). This manifests the periodicity of the potential and the infinite de- 



generacy of the vacuum. The minima of the potential ( |3.57|) occur at 
(rij + 1/2)71 for integer rij thus the degenerate vacua are specified by a set 
integers (ni, n2, n^). The rank r of B'^B is one for SU{2)/U{1) and 
SU{3)/U{2) and two for SU{A)/S{U{2) x U{2)). Thus solitons for the 
SU{4:)/ S{U{2) X f/(2)) case, interpolating between two vacua {nia,n2a) and 
{nn),n2b), are labeled by two soliton numbers Ani = rin, — riia and An2 = 
— n2a- As a specific example, we may write the B matrix of SU{3)/U{2) 

by 

B = {-(f)smT]e-'^ - <pcosr]e-"^) (3.58) 
so that the m matrix (|3.54|) is given by 



m 



I COS0 — sin sin r/e"*^ — sin^cosr^e" 

sin^sinr^e*^ cos^ ?7 + cos sin^ ?7 — cos77sin?7e*^~*"(l — cos< 

y sin0cos?7e*" — cosr7sinr7e*°~*^(l — cos0) sin^ r/ + cos cos^ r/ 

(3.59) 

Then, the potential term becomes 

Tr(^T^"^f ) = Tr(mTm-^T) = 6-9 sin^ (3.60) 
which agrees with ( |3.56|) . 



4 Dressing and solitons 

In order to find exact solutions of the SIT equations with inhomogeneous 
broadening, we first note that the linear equation ( 3.51|) admits an application 



of the dressing method. The dressing method is a systematic way to obtain 
nontrivial solutions from a trivial one. In our case, we take a trivial solution 

by 

(/ = 1 and ^ = ^° = exp[-(A - + ~\T)z - (j^-j)^ A- (4-1) 

Let r be a closed contour or a contour extending to infinity on the complex 
plane of the parameter A and G[X) be a matrix function on F. Consider the 
Riemann problem of \E'°G(A)(\E'°)~^ on F which consists of the factorization 

*°G'(A)(^°)-i = (<l>_)-i$+ (4.2) 
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where the matrix function ^, A) is analytic with n simple poles fii, Hn 

inside F and ^-{z, z, A) analytic with n simple zeros Ai, A„ outside F. We 
assume that none of these poles and zeros lie on the contour F and the 
factorization is analytically continued to the region where A 7^ /ij, Aj ; i = 
l,...,n. We normalize $+,<l>^ by = $-lx=oo ~ E)ifferentiating 

( |4.2|) with respect to z and z, one can easily show that 

B^.K^-{^) ^ (4.3) 

Since $+($_) is analytic inside (outside) F, we find that the matrix functions 
tJ and V, defined by 

U = -d^<^-^ + ^{A - + XT)^-^ - ~XT 

V = -{X - Od"^^'^ + "^Tf^-^ (4.4) 

where $ = $+ or $_ depending on the region, become independent of A. 
Then, \E' = $\E'° satisfies the linear equation; 

{d + U + XT)^ = , {d+ ^^-^^)^ = . (4.5) 

Since U, V are independent of A, we may fix A by taking A = ^. Define g 
hy g = -ff$~^|x=5 where H is an arbitrary constant matrix which commutes 
with T, T and A. Then, U and V become 

U = g-'dg + g-'Ag-iT (4.6) 
V = g-'Tg. (4.7) 

If we further require the constraint condition ( p.50| ) on such that 

{-d^^-^ + ^A^-^)h- A = 0, (4.8) 

we obtain a nontrivial solution g and \E' from a trivial one. The nontriv- 
ial solution in general describes n-solitons coupled with radiation mode. If 
G{X) = 1 in ( [4.2|) , we obtain exact n-soliton solutions. This formal procedure 
may be carried out explicitly for each cases of SIT in Sec. 3 and a closed 
form of n-soliton solutions can be obtained. In the following, we give explicit 
expressions of 1-soliton solution of two different cases and discuss about their 
physical properties. 
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4.1 Nondegenerate two- level case 

Here, we set /5 = 1 without loss of generality. The 1-soliton solution may 
be obtained either by using the above dressing method or by applying the 
Backlund transformation directly to the trivial vacuum solution [0. It is 
given by 

cos 09 = , sech(2fe — 2bCz) 



ri = {a-^)z+{a-^)Cz 

9 = - tan"^ [^^coth(262 - 2bCz)] - 2^z + 2Dz (4.9) 
where a, b are arbitrary constants and 

^ = {(a-ey + b^) ' ^ = " {^a-^'y + b^) ~ {(a-^'y + b^) ' 

(4.10) 

In terms of E as defined in (|2.13| ), 1-soliton is given by 

E = -2ib sech(262 - 2bCz)e-^'^'''-^~'^^''-^'^^~^\ (4.11) 

In the sharp line limit of the frequency distribution /(^') = ^('C'— 0? 1-soliton 
retains the same form except for the change of constants C and D, 

The name "n-soliton" is rather ambiguous since even in the sharp line limit 
the number n does not necessarily mean the topological soliton number. For 
example, when a = ^ in the above 1-soliton solution, the solution describes a 
localized pulse configuration which interpolates between two different vacua 
of (|3.29|) such that 

4>{x = -oo) = {n + i)7r, 0(x = oo) = (n + i - |^(-l)")vr. (4.13) 

Thus it carries a topological number An = (— 1)"^^6/|6| and becomes a 
topological soliton. However, when a ^, the solution interpolates between 
the same vacuum since the peak of the localized solution does not reach 
to the point where cosif = 1. Therefore, its topological number is zero. 
Nevertheless, it shares many important properties, e.g. localization, scatter- 
ing behavior etc., with the topological soliton so as to deserve the name, a 
"nontopological soliton". One of the important feature of a nontopological 
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soliton is that instead of a topological charge, it carries a [/(I) charge which 
accounts for stability of a nontopological soliton. The issue of U{1) charge 
and stability will be addressed in the following sections. It is remarkable 
that the nontopological soliton can be obtained from the topological one by 
the local vector transform ( |3.15| ). Note that even though the action ( p.l3| ) 
is invariant under the gauge transformation, the physical quatities like E, P 
and D are not invariant. The gauge choice ( |3.25| ) refers to the amount of 
detuning with a detuning parameter Au = ^. Therefore, by choosing a dif- 
ferent gauge condition through the vector gauge transformation, we obtain 
a system with a different amount of detuning. In particular, we can obtain 
a solution g{Au! = ^) of SIT theory with detuning parameter ^ from the 
solution with zero detuning through the transform 

g{Auj = = h~^g{Auj = 0)h, h = exp{^zT). (4.14) 

Or, in case of soliton solution, the topological soliton { a — C, =0) can be 
mapped to the nontopological one (a — ^ 7^ 0). For this reason, in this paper 
we will simply call n-solitons for both the topological and the nontopological 
solitons unless otherwise stated explicitly. 



4.2 Degenerate three-level case 

We first restrict the degenerate three-level system to the resonant case (Ai = 
A2 = 0, V = 0) without external magnetic field and inhomogeneous broad- 
ening. This is equivalent to the case where A = A = in (|3.20|) with 



identifications ( p.46| ) in terms of a 4 x 4 matrix g. By applying the Backlund 
transformation for the symmetric space sine-Gordon theory in Ref. |T^, we 
obtain the 1-soliton solution with the parametrization as in ( |3.54| ), 

2 

B = -2Bq tan"^ exp(26;z ^ ^ const.) = (I)Bq (4.15) 
where 6 is a constant and Bq is a constant 2x2 matrix satisfying 

BoBIBo = Bo. (4.16) 
If the matrix Bq is degenerate, i.e. det-Bo = 0, Bq can be given in general by 

with complex constants 0,61,62 satisfying |6'ip + |^2p = 1- The eigenvalues 
of BqBq are then zero and one. Therefore, up to a global SU{2) similarity 
transform of BqBq, this solution corresponds to the (1,0) or (O,l)-soliton as 
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explained in Sec. 3.3. This solution has been known as a simulton in ear- 
lier literatures. Note that the effective potential ( ^.101 ) is invariant under 
the similarity transform. This sounds to be contradictory since the similar- 
ity transform mixes two different emission modes of radiation with different 
resonance frequencies, Wha 7^ Whc, consequently different transition energies. 
However, in considering the degenerate three-level case, we have assumed 
that the oscillator strengths are equal, 

ki\dba\' = h\dtf, (4.18) 

for incident wave numbers fci, k2 and the reduced dipole momenta dba, dbc [0]- 
This assumption, together with a mediatory role of off-diagonal components 
of the density matrix, the similarity transform in fact becomes an invariance 
of the potential. 

For the nondegenerate Bq, we can take Bq as an arbitrary U (2) matrix so 
that BqBq = l2x2 and the corresponding solution is the (l,l)-soliton in Sec. 
3.3. This is different from (1,0) or (O,l)-soliton and can not be reached by 
the similarity transform since the similarity transform preserves eigenvalues 
of BqB^. In fact, they are topologically distinct and separated by an infinite 
potential energy barrier. Finally, physical quantities can be obtained from g 
through the identification ( |3.2CI|) . Explicitly, we find E, P and M in ( ^.46| ) 
to be 

2_ 

E = iBodcj) = —2ir]Bosech.(2riz H — z + const.) 

V 

P = -25osin20, 

M = -A^ = -212x2 cos 20 (4.19) 

respectively. Inclusion of detuning and external magnetic effects can be done 
easily by a gauge transform; 

E H^^EH2, M H:[^MHi, P H:[^PH2, N H^^NHi (4.20) 

where Hi, H2 are given by = H^^dHi, A2 = H^^dH2 for A^, A2 in (|^). 



5 Symmetries 



The effective field theory action ( p.l3|) reveals various types of symmetries of 



SIT. As an integrable field theory, at least classically, it possesses infinitely 
many conserved local integrals. These conservation laws can be extended 
to the inhomogeneously broadened case without difficulty. Some explicit 
local integrals for the nondegenerate two-level SIT have been given in earlier 



literatures [^. Here, we present a systematic way to find local integrals 
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for the general G/H case and present a few explicit examples. Another 
important symmetry of ( ^.1^ ) is the global U{1) axial symmetry; g — * hgh 
for a constant element h G U{1) C H. This type of symmetry has been 
previously unknown and is one of the outcomes of our effective field theory 
approach. We show that the U{1) symmetry and its corresponding charge 
plays a crucial role in pulse stability problem. 

Besides the aforementioned continuous symmetries, the action ( p.l3[ ) also 
possesses two distinct types of discrete symmetries; the chiral and the dual 
symmetries. Discrete symmetries relate two different solutions of the SIT 
equation. In particular, the dual symmetry relates the "bright" soliton with 
the "dark" soliton of SIT. These discrete symmetries reflect the ubiquitous 
nature of the action (|3.13|) without the potential term, as an action for the 
coset conformal field theory and also a bosonized version of 1+1- dimensional 
free fermion field theories. 



5.1 Conserved local integrals 

In the previous section, we have shown that the linear equation with a spec- 
tral parameter yields exact soliton solutions through the dressing procedure. 
The same linear equation can be employed to construct infinitely many con- 
served local integrals. In this section, using the linear equation as well as 
the properties of Hermitian symmetric spaces, we present a systematic way 
of finding such integrals for various cases of SIT introduced earlier. 

We first recall some relevant mathematical facts on Hermitian symmetric 
space A symmetric space G/H is a coset space with the Lie algebra 



commutation relations among generators of associated Lie algebras such that 

[h , h] C h , [h , m] C m , [m , m] C h , (5.1) 

where g and h are Lie algebras of G and H and m is the vector space 
complement of h in g, i.e. 

g = h©m. (5.2) 

Hermitian symmetric space is a symmetric space equipped with a complex 
structure. In general, such a complex structure is given by the adjoint action 
of To on m up to a scaling, where Tq is an element belonging to the Cartan 
subalgebra of g whose stability subgroup is H. In our case, Tq is precisely 
the T-matrix given in Sec. 3. Namely, with a suitable normalization of T, 
we have 

T G h , [T , h] = and [T , [T , a]] = -a for any a G m. (5.3) 
We decompose an algebra element ^ § according to ( |5.2| ), 

i' = ^h+ ^m- (5.4) 
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Such a decomposition could be extended to a group element \I' G G = SU (n) 
if we substitute the commutator by a direct matrix multiplication and add 
an identity element Hq = I to the subalgebra h, i.e. 

^ = ^h + C ^h, ^/^^m C ^ m C ^ft. (5.5) 

In other words, any unitary n x n matrix can be expressed as a linear com- 
bination of SU (n) generators and the identity element /iq such that 

dim h dim m 

^ = ^^ + ^,„= J2 E (5-6) 

a=0 6=1 

In order to solve the linear equation ( |3.51| ) recursively, we expand \E' in terms 
of a power series in A, 

oo 1 dim h dim m 

^exp{-\Tz)=Y,j:<i>i, where $i = ^ E ^■mb = Vi + Qi. 

(5.7) 

We also introduce the abbreviation: 

S = U = g'^dg + g'^Ag - C m 
(V)^ ^ {g-'Tg)^ = {g-'Tg{-^'y) = Di + Pi, A C h, C m(5.8) 

so that the linear equations become 

{d + S)^,-[T , $.+i] = (5.9) 

and 

i-l 

+ + P^-l-l)^l = 0. (5.10) 



1=0 

Then the m-component of (|5.9|) is 

dQi^i + SVi-i - [T , Qi] = 0, (5.11) 

which can be solved for Q,; by applying the adjoint action of T, 

Qi = -[T , 9Q,_i] - [T , S]V,^i. (5.12) 

Vi can be solved similarly from the h-component of (|5.9|) and (|5.10| ) such 
that 

i—l 

V, = - [ SQdz -J2 [ (A-z-iP; + P^-l-lQl)dz. (5.13) 

1=0 
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These recursive solutions of Vi and Qi can be determined completely with 
appropriate initial conditions. For example, if we choose an initial condition 
which is consistent with the recursion relation for i < 0, 

V, = I , Qo = 0, (5.14) 

we find for the first few explicit cases in the series, 

Vi= j £[T , £]dz - I Dodz, Qi = -[T , £] (5.15) 

and 

P2 = / (^5^ + ^[T , £]Vi)dz + I (-Di - DoVi + Po[T , £])dz 
Q2 = -d£ - [T , £]Vi. (5.16) 

Finally, the consistency condition: ddVi = ddVi gives rise to infinitely many 
conserved local currents, 

J, = dV. = -£Q,, J, = dV^ = - ^(A-^-iP/ + Pi^i-iQi), (5.17) 

1=0 

satisfying dJi = dJi. A few examples are 

Ji = £[T , £] , Ji = -Do (5.18) 

J2 = £d£ + £[T , £]Vi , J2 = -Di - DoVi + Po[T , £]. (5.19) 

The first current Ji, Ji gives rise to the energy conservation law. With the 
repetitive use of the properties of the Hermitian symmetric space, it can 
be easily checked that these conservation laws are indeed consistent with 
the equations of motion ( |3.52D ( ^^53|) , which in the present convention take a 
particularly simple form: 

d£-[T , Po] = 
dD, + [£ , P,] = 
dP, + [£ , A] - [T , P.+i] = 0. (5.20) 

In general, the conserved current contains nonlocal terms. These nonlocal 
terms may be dropped out by taking the T-component of the currents. For 
instance, the T-component of the "spin-2" current conservation is 

dTr{T£d£) = dTT{TPo[T , £] - TD{) (5.21) 
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which obviously does not contain nonlocal terms. In order to demonstrate 
the above procedure more explicitly, we take the SU{2)/U{1) case as an 
example. In this case, relevant 2x2 matrices are given by 

^-(l -n^l. T^U: 1) (8.22) 



g-'Tg) = (g-'Tg(-n') ^ -i . , , , , , 



,E* y ' 2 V -i 

and we introduce the notation 

Let 

vl/exp(-AT;.) = ^^; ^ ^0 (5.24) 

so that the linear equation changes into 

(^+(^* "o^))*^-[^' *m] = (5.25) 

and ^ 

d<^>, + J2{9''fg) <^i^0. (5.26) 

1=0 * 

These equations can be solved iteratively in component, 

Qi = hdqi.i-Esi.i) (5.27) 
Zl 

Ti = -^(9r,_i + £;*p,_i) (5.28) 



Pi ^ f Endz + iY.f {Di_i_ipi + Pi_i_iri)dz (5.29) 

= - [ E*qidz + iY, [{-Di-i-iSi + P*_i_,qi)dz (5.30) 

together with the initial conditions: 

Po = So = -2i, ro = gro = 0. (5.31) 

The consistency: ddpm — ddp^ — leads to the infinite current conservation 
laws, dJi + dJi — 0, or 

id ^{Di_i_,pi + Pi^i^iri) - d{En) = 0. (5.32) 

1=0 



31 



The consistency condition, ddsi — ddsi = 0, gives rise to the complex conju- 
gate pair of ( p. 321 ). A few exphcit examples of conserved currents are 

Ji = -2Do 

Ji = EE* (5.33) 

Js = AiDi-2PoE* 

J2 = EdE* (5.34) 

J3 = -2PodE* + 8D2 + UE*Pi 

J3 = Ed^E* + {EE*f (5.35) 

J4 = -IGiD^ + 8E*P2 + AiPidE* - 2Pod'^E* - 2PoE*\E\'^ 

J4 = Ed^E* + \E\'^d\E\'^ + 2E\E\'^dE* (5.36) 

5.2 ^(1) symmetry 

One of the properties of the Wess-Zumino-Witten action ( p.9|) is the global 
axial vector gauge symmetry, i.e. 

Swzwifgf) = Swzwig) (5.37) 

for a constant f E G. The extra terms in the deformed action ( |3.13| ) break 
this global axial symmetry in general. Nevertheless, there remains at least an 
unbroken [/(l)-axial symmetry given hj g ^ hgh for h = exp(7T) G U{1). 
This t/(l)-invariance results in a conserved charge according to the Noether 
method. Even though a general expression for the conserved charge should 
be possible, in practice it requires an explicit parametrization of the group 
variable g. Thus, for the sake of brevity, here we will restrict ourselves only 
to the SU{2)/U{1) case. In this case, the global axial transform is given by 

7] ^ 1] + '-f for 7 constant. (5.38) 

The corresponding Noether currents and the associated axial charge are 

sin Lf sin Lf 

/oo _ 
dx{J + J) (5.39) 
-00 

where, owing to the complex sine-Gordon equation (|2.14|) - (|2.16|) , J and J 
can be shown to satisfy the conservation law: dJ + dJ = 0. In particular, 
the axial charge of the nontopological 1-soliton in ( [4.9|) is 

Ql-sol = ctan-i^^. (5.40) 
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The charge of the topological soliton is not well defined ||T9|. Stability of 
nontopological solitons can be proved either by using conservation laws in 
terms of charge and energy as given in [|1^, or by studying the behavior 
against small fiuctuations which will be explained in Sec. 6 p2| . 

5.3 Discrete symmetries 



Besides continuous symmetries, the action ( p. 131) also reveals discrete symme 



tries of SIT, the chiral symmetry and the dual symmetry. They are manifested 
most easily in the gauge where A = A = Q. Extensions to different gauges, 
e.g. the off- resonant case which requires a different gauge fixing as in ( |3.25| ), 
can be made by the vector gauge transform in ( p.l5|) . 



One peculiar property of the action ( p.l3| ) is its asymmetry under the 
change of parity z ^ z. This is because the Wess-Zumino-Witten action 



(|3.9|) is a sum of the parity even kinetic term and the parity odd Wess- 
Zumino term thereby breaking parity invariance. In the SIT context, broken 
parity is due to the slowly varying enveloping approximation which breaks 
the apparent parity invariance of the Maxwell-Bloch equation. Nevertheless, 
the action (|3.13|) is invariant under the chiral transform 

z ^ z , g ^ g^^ ( 01 r] ^ —r] , y9 ^ —ip) (5-41) 

which may be compared with the CP invariance in the context of particle 
physics. Thus, parity invariance is in fact not lost but appears in a different 
guise, namely the chiral invariance. This chiral symmetry relates two distinct 
solutions, or itgenerates a new solution from a known one. For example, 
under the chiral transform ( |5.41| ), the 1-soliton solution ( [4.9| ) in the resonant 



case = 0) becomes again a soliton but with the replacement of constants 
a, h by 

-^-^-'"-^r (5-42) 

This implies the change of pulse shape and the change of pulse velocity by 
f ^ c — f . The current and the charge also change into 

J ^ - J , J ^ - J , -Q. (5.43) 

It is remarkable that the velocity changes from v to c — v unlike the usual 
parity change f — > — f . 

The other type of discrete symmetry of the action (p.l3|) is the dual 
symmetry of the Krammers-Wannier type: ||T9l 



/3 ^ -/3 , g^iag (5.44) 
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where a is a constant matrix with a property aT + Ta = 0. For example, 
0" = (Ti of Pauh matrices for the SU{2)/U{1) case. This rather unconven- 
tional symmetry, also the name, stems from the ubiquitous nature of the 
action (|3.13| ), i.e. it also arises as a large level limit of parafermions in statis- 



tical physics and the above transform is an interchange between the spin and 
the dual spin variables [0 . In general, the change of the sign of (3 makes the 
potential upside down so that the degenerate vacua becomes maxima of the 
potential and vice versa. Therefore, the dual transformed solutions are no 
longer stable solutions. This allows us to find a localized solution which ap- 
proches to the maximum of the potential asymptotically (so called a "dark" 
soliton). In practice, the dark soliton for positive /3 can be obtained by re- 
placing P — > —P , z —z in the "bright" soliton (usual one) of the negative 
P case. For example, we obtain the dark 1-soliton for the SU{2)/U{1) case 
as follows: 

cosipe"^'" = ^ tanh{2bz + 2bCz) ' 



a-O^ + b^ ^{a-if + b^ 

2{a-i){z-Cz)-2iz. (5.45) 



6 Stability 

The physical relevance of a soliton number is that it accounts for the stabil- 
ity of solitons against "topological" (soliton number changing) fluctuations. 
Note that any finite energy solution must approach to one of the degenerate 
vacua at both ends of the thus carries a specific soliton number. In 

general, solutions are made of soliton, antisoliton, breather and radiation 
modes where each modes carry soliton number 1, -1, and respectively. 
The overall soliton number is a conserved topological charge which satisfies 
the so-called superselection rule, i.e. it cannot be changed during any phys- 
ical processes due to the infinite potential energy barrier between any two 
finite energy solutions with different topological charges. This infinite en- 
ergy barrier results from the infinite length of the x-axis despite the finite 
potential energy density per unit length. In the case of sine-Gordon solitons, 
the "area" of the McCall and Hahn's theorem just enumerates the soliton 
number, i.e. 

/OO /"OO 
2Edt = / 2dipdt = 2(p{t = oo) - 2(p{t = -oo) 
-OO J —oo 

= 2<p{x = -oo) -2<p{x = oo){= -2TTAn) (6.1) 

where the equality in the last line holds for localized solutions. However, this 
coincidence is only for a specific case where E is real and inhomogeneous 
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broadening is ignored. As was shown in Sec. 2, the inclusion of a phase 
degree of freedom into E makes E complex so that the time area of 2E 
becomes complex too, thus it loses its meaning as a soliton number. 

Inhomogeneous broadening also spoils the notion of a soliton number. In 
Sec. 3, we have noted that the effective potential variable g becomes a micro- 
scopic variable depending on the frequency ^ and atomic variables P and D 
also depend on ^. This implies that the soliton number is also a ^-dependent 
microscopic variable. However, the macroscopic field E, also given in terms of 
g, is independent of ^. In other words, the microscopic variable g arranges the 
^-dependence in such a way that the resulting E field becomes independent 
of ^. This is exemplified by the 1-soliton solution in ( [4.9D -( ^?TTD . Thus, inho- 
mogeneous broadening in general requires to be a function of "frequency 
^ averaged" coefficients so that it does not carry a topological soliton num- 
ber. In this regard, the McCall and Hahn's area theorem is different from 
the topological stability argument and 2mi pulses should be distinguished 
from n-solitons. It is remarkable that the area theorem provides a stability 
argument even in the absence of the topological argument. In fact, the proof 
of the area theorem relied crucially on the averaging over the frequency ^ 
of detuning associated with inhomogeneous broadening. However, one seri- 
ous drawback of the area theorem is that it applies only to a very restricted 
case, i.e. the case of real E (neglecting frequency modulation effect) and the 
symmetric frequency distribution. Presently, a more general area theorem 
including frequency modulation is not known and finding such a theorem 
would be one of the most important problem in the theory of self-induced 
transparency. 

For the rest of the section, we consider a restricted case which however 
generalizes the area theorem to a certain extent. We consider the SU{2)/U{1) 
case without inhomogeneous broadening and prove the area stability around 
the 1-soliton solution (topological or nontopological) using a perturbative 
argument similar to the one given in pO|. This shows that the reshaping of 
a pulse also occurs in the off resonant case (a — ^ 7^ 0) but slower than in the 
resonant case. We argue that the U{1) charge of optical pulses introduced 
in earlier sections could be used also in generalizing the area theorem. In 
particular, we prove the U{1) charge stability using the local conservation 
law. These results are in good agreement with numerical results obtained 
earlier 123 . 



Consider the bright soliton (|4.9| ) which describes a coherent pulse propa- 
gating in the attenuator. In order to overcome the difficulty with the com- 
plexity of the usual area function for complex E, we regard of the complex 
sine-Gordon equation as a "modified" area function. Without loss of gener- 
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ality, we may take the asymptotic time behavior of if by 

= — oo, x) = — 7r/2, ip{t = oo,x) = Ti /2 for a = ^ 
ip{t = —oo,x) = ip(t = oo,x) = —tt/2, for 07^^, (6.2) 

so that the modified area _ 

A= / 2d<^dt (6.3) 



of the topological soliton (a — ^ = 0) is 2tt while that of the nontopological 
soliton is zero. Now assume that the system was initially in the vacuum state 
{(f(t = —oo,x) = —7t/2). If the solution is perturbed around a soliton so 
that near the trailing edge of the pulse (t » 1), the modified area function 
changes hj Sip = e for small e, i.e. (p{t » 1) = ±7r/2 + e. Then, the complex 
sine- Gordon equation for t >> 1 becomes 

where we have neglected the variation of rj which contributes only to the 
order e^. This shows that if the modified area is greater than 27r (or zero) 
by the amount e > 0, then Bdip < so that the field dip at the trailing edge 
tends to decrease along the z = x/c axis to recover a total modified area 
271 (or zero). On the other hand, if the perturbation is such that e < 0, 
then ddip > and the field at the trailing edge increases. Therefore, the 
total modified area tends to remain 2tt or zero. Moreover, ( |6.4|) shows that 
the recovery of area is faster in the resonant case (a = ^ ) than in the off- 
resonant case (a 7^ ^) which agrees precisely with the numerical result 



In fact, the recovery of the modified area is acompanied by a stronger pulse 
reshaping recovery to that of a soliton. Here, we only point out that the 
stability of a soliton shape against fluctuations which preserve the modified 
area may be shown by making a modification of the Lamb's proof using 
the Liapunov function pOf , and also by proving the stability of higher order 



conserved charges as discussed below. In case of the dark soliton ( p.45| ) 
which describes a coherent pulse in the amplifier, the system is initially in 
the upper level {ip(t = —oo,x) = 0). In such a case, one can easily see that 
the above consideration of perturbation results in the instability of a soliton 
configuration. 

It was suggested previously that pulse reshaping arises due to the stability 
of all higher order conserved integrals against small fluctuations and each 
integrals behaves like an "area" |Q. In order to explain the stability of 
conserved integrals, we apply the method introduced in to the f/(l)-charge 
conservation law: 

d d 

— [cot^ ^{d + d)ri]+ c— [cot^ (pdr]] = 0. (6.5) 
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Introduce the U{1) charge in terms of a "time area", 



/ + 00 
[ccotV^^]^^^, (6.6) 
-oo 



so that fB.SI) becomes 



dA 

dx 



- cot' V{d + d)v\lZtZ (6.7) 



The boundary contribution from the 1-sohton is zero thus the ?7(l)-charge 
A{x) is conserved in space, i.e. dA/dx = 0. The physical meaning of the U{1) 
charge is clear; since dr] = a — ^ = a + Wo — w for the 1-soliton which expresses 
frequency detuning and cot^y? is peaked around the soliton, the [/(l)-charge 
measures precisely the amount of detuning. We caution that, as already 
noted by Lamb , the carrier frequency Wo shifts to Wo + a in the presense 
of a soliton so that the detuning parameter is given hj a + Wo — w. Now, if 
the solution is perturbed around the soliton such that near the trailing edge 
of the pulse, 

(l){t»l,x) = ±7r/2 + e(x) 

Vit»l,x) = (a-0{t-^-j^-^^^^)+6{x) (6.8) 

for small parametric functions e{x) and 6{x). To the leading order, the 
variation of the f/(l)-charge then becomes 

d6A , , ( 1 \ o , , 

This shows that the detuning by a higher frequency, i.e. a — ^ > reduces A 
for increasing x while the lower frequency detuning does exactly the opposite. 
Since the conserved charge A of the 1-soliton is ctan~^[|&|/(a — 0]' can 
be concluded that the absolute value of A{x) decreases monotonically along 
the X-axis and converges to a constant charge of the soliton. Note that the 
recovery of |y4| value to that of the soliton is slower than the area case since 
it is of the order e^. The decreasing behavior of |y4| is in good agreement with 



the numerical work which showed that the frequency of the optical pulse 
is pulled towards the transition frequency and reaches to a constant value 
along the x-axis. Thus, the ?7(l)-charge stability provides a generalization of 
the area theorem in the presense of frequency detuning. This type of stability 
argument might be extended to the case of other conserved integrals. 

In the presence of inhomogeneous broadening, the U{1) conservation 
breaks down unlike the local conservation laws in Sec. 5.1. It introduces an 
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anomaly term M in the ?7(l)-current conservation such that dJ + dJ = M 
for J, J in ( |5.39| ) and 



M = 2 cot v9[ cos(6' - 2t]) < sm{9 - 2r/) sin > 

- sm{9 - 27]) < cos{9 - 2r/) sm2^ > -(cot^ ^dr] + ^d9)d^ ](.6.10) 



This anomaly vanishes in the sharp line limit due to the constraint ( |2.16| ). It 
also vanishes in the case of 1-soliton and the charge remains conserved. This 
behavior may be compared with the conserved area of topological solitons 
in the presence of inhomogenous broadening. The area theorem of McCall 
and Hahn proves that inhomogeneous broadening changes the pulse area 
until it reaches to those of 2n7i pulses. It remains as an important open 
question whether one could prove a generalized area theorem of pulse stability 
including frequency modulation by making use of U{1) charge and anomaly. 



7 Discussion 

In this paper, we have presented a field theoretic formulation of SIT. Various 
cases of SIT have been associated with specific symmetric spaces G/H and a 
systematic group theory approaches to the SIT problem has been made. In 
doing so, the introduction of a matrix potential variable g was an essential 
step. One immediate question is about the generality of such a group vari- 
able in the description of nonlinear optics problems. In the nondegenerate 
two-level case, the SU (2) variable g was introduced by solving the constraint 
which expresses the conservation of probability. In more general multi-level 
cases, the complexity of atomic states requires bigger groups G than SU{2). 
However, not all degrees of freedom associated with bigger groups become 
physical degrees of freedom, they are constrained by the Bloch equation which 
requires a restriction of g e.g. as in ( |3.26|) . That is, the microscopic atomic 
variables couple to the macroscopic electric field selectively through dipole 
transitions which suppresses the excitation of certain degrees of freedom. 
This reduction of variables has been accomplished by considering the G/H 
coset structure. It should be emphasized that the effective field theory we 
introduced is not an ordinary sigma model where g takes value in the coset 
G/H. In the present coset structure is imposed only through a non- 

trivial Lagrange multiplier.^ Unfortunately, this coset theoretic description 
is not valid for all cases of SIT. For example, the particular way electric field 
interacts with two-level degenerate states in the 2^2 transition does not 
allow a nice coset structure, but only described by the double sine-Gordon 

^ More correctly, it corresponds to the coset LG/LH of loop groups LG and LH rather 
than G/H, see Ref. || 
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equation when further restrictions are made. Thus, a systematic apphcation 
of group theory to the general SIT problems still remains as an interesting 
open question. 

On the other hand, the group theoretical approach is not restricted to 
the SIT problem only. The nonlinear Schrodinger equation, which is the 
governing equation for optical soliton communication systems, can be also 



generalized according to each Hermitian symmetric spaces In fact, both 
SIT and the nonlinear Schrodinger equation share the same Hamiltonian 
structure and can be combined together. We will consider this case and its 
physical applications in a separate paper. 

Finally, our approach provides a vantage point to the quantum SIT prob- 
lem as well as the quantum optics itself. A direct quantization of SIT using 
the quantum inverse scattering has been made by Rupasov and a localized 
multiparticle state has been found and compared with a quantum soliton 

. Our group theory formulation of SIT suggests an alternative, yet more 
systematic way of quantizing the SIT theory according to each specific coset 
structures. More generally, our approach can be extended to other quantum 
optical systems. Thes works are in progress and will be reported elsewhere. 

ACKNOWLEDGEMENT 

This work was supported in part by the program of Basic Science Re- 
search, Ministry of Education BSRI-96-2442, and by Korea Science and En- 
gineering Foundation through CTP/SNU. 

References 

[1] S. L. McCall and E. L. Hahn, Phys. Rev. Lett. 18, 908 (1967); Phys. 
Rev. 183, 457 (1969). 

[2] A. I. Maimistov, A. M. Basharov, S. O. Elyutin and Yu. M. Sklyarov, 
Phys. Rep. 191, No. 1 , 1-108 (1990). 

[3] G. L. Lamb, Jr., Phys. Rev. Lett. 31, 196 (1973). 

[4] M. J. Ablowitz, D. J. Kaup and A. C. Newell, J. Math. Phys. 15, 1852 
(1974). 

[5] A.I. Maimistov, Opt. Spektrosk. 57, 564 (1984). 

[6] A.M. Bashrov and A.M.Maimstov, Sov. Phys. JETP 60, 913 (1984). 

[7] A.M. Bashrov and A.M.Maimstov, Sov. Phys. JETP 67, 2426 (1988). 



39 



Q-H. Park and H. J. Shin, Complex sine-Gordon theory for coherent 
optical pulse propagation, [hep-th/950912U . 



[9 
[10 

[11 
[12 

[13: 

[14 
[15 
[16 

[17; 

[18 
[19 
[20 
[21 

[22 
[23 
[24 
[25 



F. Lund and T. Regge, Phys. Rev. D14, 1524 (1976). 
K. Pohlmeyer, Commun. Math. Phys. 46, 207 (1976). 

F. Bloch, Phys. Rev. 70, 460 (1946). 

E. Witten, Commun. Math. Phys. 92, 455 (1984). 

D.Karabali, Q-H. Park, H. J.Schnitzer and Z.Yang, Phys. Lett. B216, 
307 (1989); K.Gawedski and A.Kupiainen,, Nucl. Phys. B320, 625 
(1989). 

I. Bakas, Int. J. Mod. Phys. A9, 3443 (1994). 
Q-H. Park, Phys. Lett. B328, 329 (1994). 

Q-H. Park and H. J. Shin, Phys. Lett. B347, 73 (1995); Nucl. Phys. 
B458, 327 (1996); T. J. Hollowood, J. L. Miramontes and Q-H. Park, 
Nucl. Phys. B445, 451 (1995). 

I. Bakas, Q-Han Park and H. J. Shin, Phys. Lett. B372, 45 (1996). 

A.M. Basharov, Sov. Phys. JETP 70, 94 (1990). 

Q-H. Park and H. J. Shin, Phys. Lett. B359, 125 (1995). 

G. L. Lamb, Jr., Rev. Mod. Phys. 43, 99 (1971). 

S. Helgason, Differential geometry. Lie groups and symmetric spaces. 
2nd ed. New York; Academic Press 1978. 

Q-H. Park and H. J. Shin, to appear. 

J.C. Diels, Phys. Lett. A 31, 111 (1970). 

Q-H. Park, Phys. Lett. B223, 175 (1989). 

A.P. Fordy and P.P. Kuhsh, Commum. Math. Phys. 89, 427 (1983); P. 
Oh and Q-H. Park, More on generalized Heisenberg ferromagnet models, 
|hep-th/960403Ti 

V.I. Rupasov, Pis'ma Zh. Eksp. Teor. Fiz. 36, 115 (1982); Sov. Phys. 
JETP 56, 989 (1982). 



40 



